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ABSTRACT 

A  nathemtloal  model  of  a  general  marine  gae  turbine  ayetem  vaa 
dereloped  in  effort  to  provide  an  analytical  representation  from  which 
the  dynamics  equations  of  specific  systems  are  dedueible.  The  model  ia 
eoigpriaed  of  a  pair  of  coupled,  non-linear,  first-order,  ordinary  differen¬ 
tial  equations.  Supplemental  data  is  requiird  to  effect  speelallBation  of 
these  equations  to  establish  representations  of  specific  systeniB,  These 
data  are  derivable  from  the  steady-state  characteristics  of  the  engine 
and  the  dynamic  responses  of  the  fuel  control  systom  to  select  input 
excitations. 
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SUIMARY  PAOE 


at  ..FEg^MI 

To  dflwiop  k  aaithemtlcal  rRpreROntntlon  or  a  general  marine  gaa 
turbine  eysteia. 

Fljdlngfl 

The  equations  which  govern  the  transient  behavior  of  such  a  system  can 
be  derived  from  phyeical  principles  independent  of  specific  system  design 
information.  Consequently,  these  equations  can  be  applied  to  a  variety  of 
gae  turbine  engines  for  purposes  of  simulation  or  control  system  develojinent. 
ReeOBwapdatloafl 

Specialise  these  equations  to  obtain  a  modal  of  the  LM^2^00  engine. 
Develop  e  conputar  sloiulatlon  of  the  model  system  in  preparation  for 
correlation  studies  with  the  I'eal  system. 
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administrative:  intormation 

Authority  for  the  research  deocrllx^d  in  this  report  is  contained  in 
NAVSKCPHIIJiDIV  Semi-Annual  Program  Suiuniery  of  1  November  1972,  The 
RDT&E  Subpro.Jeot  No.  assigned  was  S-F 32-10!> ,  Task  12f)Q4;  the 
MAVSECPHILADIV  Project  No.  is 

This  is  a  final  report,  the  seco* d  of  a  series^^^  of  reports  describing 
etadies  related  to  the  dynamics  of  a  marine  propulsion  gas  turbine  system. 
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REPORT  OF  INVESTIGATION 

lirrRODUCTION 

The  need  for  a  generalized  analytical  model  of  a  gae  turbine  aystem 
is  generated  by  the  variety  of  engines  being  used  or  proposed  for  naval 
prop^iiaion  plants.  Representations  of  specific  systems  are  special iaations 
of  this  iDO'lel  achieved  by  imposing  the  requisite  mathematical  constraints. 

The  resu-ltlng  equations  can  be  mechanized  on  a  computer  to  serve  as  a 
dynamic  test  system  for  research  being  performed  under  related  projects. 

The  DIVISION  is  currently  engaged  in  a  ..umber  o<‘  projects  concerned 
with  the  dynamic  response  characteristics  of  gas  turbine  systems. 

Design  of  a  control  system  for  a  high  .'.peed  waterbrake  has  been  undertaken 
in  an  effort  to  simulate  the  dynamic  loading  encountered  by  gas  turbines 
drlring  controllable  pitch  propellers.  Preliminary  design  evaluations  will 
be  conducted  on  a  computer,  end  therefore  will  require  a  congniter  synthesis 
of  the  plant.  This  simulation  may  also  be  used  for  developmental  testing 
of  a  method  designed  to  determine  the  dynamic  eharaeterlstlcs  of  nonlinear 
systems.  The  method  is  em  analytical  procedure  that  requires  dynaaio  rasponse 
data  acquired  by  linear  test  methods  for  the  construction  of  a  system 
representation.  The  nonlinear  characteristics  of  the  gas  turbins  syeteni 
will  require  the  application  of  the  full  machinery  of  this  analysis  to  the 
model  truusient  -uest  data  and  thereby  e^diaust  possible  problem  areas  in 
the  analysis. 
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DffiCKIPflON  OF  SYSTD4 

The  erglDB  ia  conjKaaed  of  two  major  units,  an  axial  flow  gas  generator 
and  a  free  (or  power)  turbine.  The  external  casings  of  Uie  gas  generator 
and  free  turbine  are  connected  together  to  form  a  continuous  duct  for  the 
paseage  of  gases  through  the  engine.  Coupling  between  the  gas  generator 
and  power  turbine  is  aarodynamlo.  The  gas  generator  oonsists  of  two  nulti- 
stage  congireasors ,  a  combustion  section  and  two  multi-stage  reaction  turbines. 
A  concentric  siiaft  ia  used  to  connect  the  coigpreseors  to  their  corresponding 
turbines. 

Systeo  input  energy  is  supplied  by  fuel  puap  and  controlled  jy 
a  multi-input,  single-output  fuel  controller.  The  controller  delivers 
part  of  the  fuel  to  the  engine  and  routes  the  rensiinder  back  to  the  pump 
interstage.  Rngine  acceleration  are  initiated  at  the  throttle  (one  of  the 
several  controller  inputs)  and  follow  schedules  defined  by  the  controller 
which  la  prograwaed  to  avoid  overte^perature  and  surge  limits  of  the  gas 
generator.  The  controller  also  includes  a  deceleration  schedule  t^mt 
providee  for  blowout  free  decelerations. 
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THE  GAO  TURBINE 

The  steady-state  characteristics  of  the  power  turbine  define  a  mathematical 
relationship  which  establi;ihe6  the  torque,  Qp, developed  by  the  power  turbine 
as  a  function,  f,  of  fuel  rate,  w^-.and  power  turbine  speed,  Np. 

Qp  =  f(Wf,  Np)  lu  ft  (1) 

Tne  coupling  between  gas  generator  turbines  is  aerodynamic  rather  theui 
mechanical,  however,  in  the  si^cady-state ,  the  fuel  rate  alone  prescribes 
the  speeds  of  the  compressor-turbine  pairs  bo  that  a  fixed  relationship 
exists  between  these  two  different  speeds.  Moreover,  the  unoccupied  woluoe 
between  the  high  and  low  speed  turbines  is  small.  Thus,  during  aceeleratlona, 
the  speeds  of  the  conpreasor-turbine  pairs  tend  to  change  in  phase  In 
accordance  with  their  steady-state  correspondence.  Let  and  Nj^  be  the 
respective  speeds  of  the  low  and  high  speed  turbines.  Then,  during  steady- 
state  or  transient  conditions,  the  current  speed,  Nj^  determines  the  speed 
in  the  manner  prescribed  by  the  foilowir»g  relationship. 

Hi  “  g(Nh)  rpm  (2) 

Suppose  the  gas  generator  system  resides  at  a  state  characterised  by 
the  fuel  rate,  WfQ  and  the  power  turbine  speed,  Np^.  At  the  instant, 
t  ®  0,  the  fuel  rate  is  inatantaneo'jsly  changed  to  Wf.  After 

the  time  interval,  dt,  this  change  in  input  energy  gives  rise  to  an  increase 
in  the  rotational  energy  of  the  entire  system.  If  dMi ,  dl^,  and  dM-  are 
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the  accoapAivying  speed  changes  In  the  low  epeed,  high  speed,  6uid  power 
turbine  respectively,  then  the  total  stored  energy  increase  dEg, 
achieved  in  the  interval,  dt,  is  given  by  the  following  equation: 

dEg  -  0.011  (Ji  Ni  dNi  +  +  JpNpdNp)  lb  ft  (3) 

The  constants  Ji,  Jh,  Jp  are  the  polar  nooents  of  inertia  of  the  three 
Dschanleally  independent  rotating  assemblies.  Equation  (2)  can  be 
used  to  eliminate  Mx  and  dli]^  from  aquation  (3). 

dEg  -  0,011  JpHpdNp  O.Cll  (JhNh  +  g(Hh))  dN  lb  ft  (4) 

Assume  that  gas  ganerator  speed,  N^,  and  fuel  rate  have  a  one-to-one 
correspondence,  h,  in  the  steady-state. 

Wf=h(Nh)  iVhr  (5) 


Let  the  specific  heat  release  of  the  fuel  bey.  Then  the  dimensionless 
apaeifie  fuel  consumption,  s,  is  given  by  the  equation: 


,  -  p.o';  — C 


(6) 


Eliidnate  in  equation  (6)  by  means  of  the  relationship  implied  by  equations 
(3)  and  (l)  to  obtain  s  as  a  function  of  H_  and  Njj  (  ^  ®  constant). 


.07  y  h(Nh) 

NpQp(h(llh),lip) 
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The  epeclflc  fuel  oonsump'tion  is  the  ratio  of  the  chemloal  energy 
input  to  the  mechanical  energy  output  end  aa  such  iopllee  that  storage  of 
an  additional  JpNpdNp  unite  of  rotational  energy  requires  the  release 
of  8  Jp  Np  dNp  units  of  chemical  energy.  Define  e  by  the  equation 

e  *  i-''  (8) 

r 

and  regard  e  as  a  dimensionless  variable  efficiency.  Then  the  chendcal 
energy  increase,  dEc, required  to  achieve  a  rotational  energy  change,  dig, 
is  given  by  the  following  equation: 

where  dEgp  is  the  additional  energy  stored  in  the  pcsrer  turbine  during  the 
interval,  dt.  Eliminate  dEgp  and  dEg  from  equation  (9)  and  substitute  the 
quantity,  0.2168'>'Aff  dt  for  dZc.  Divide  the  resulting  equation  by  dt  to 
determine  that; 

19.7  Wf  -  i  JpNpNp  ^  (JhNh  *  Jig'  (Nh))Hh  lb  ft  (10) 

The  tenris ,  A  Wf  la  the  difference  fuel  rate  defined  by  the  equation: 

A  Wf  (t)  -  W^{t)  -  Wf^s(t)  Ib/hr  (11) 

where  Wf(t)  is  the  actual  time-dependent  fuel  rate  and  is  the  quaal- 

etatlc  fuel  rate  determined  by  the  following  procedure.  Solve  equation 
(1)  for  Wf  and  obtain: 
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Wf.  -  C  (Op,  Np)  iVhr  (12) 

Equate  Qp  to  is  the  load  torque)  and  Wp  to  and  obtain  the  resxilt: 

Wfqe  ■  C  (Qi,  Np)  Ib/hr  (13) 

Finally,  substitute  equation  (13)  into  equation  (11)  and  replaot  A  Wf  In 
equation  (10)  with  this  result. 

19.7y  (Wf(t)  -  C  (Qi,  Np))  « 

J  JpNpNp  +  (JhNli  Jl6'(Wh))  lb  ft  (U) 

« 

The  acceleration,  Np  of  the  power  turbine  is  inversely  proportional 
to  the  total  polar  moioant  of  Inertia  (consisting  of  the  power  turbine  and 
load)  and  directly  proportional  to  the  net  torque,  A  pQ  acting  on  the 
turbine  shaft, 

Np  -2^  ^pQ  rpm/eec  (15) 

Jp 

Moreover,  the  net  torque  is  comprised  of  the  applied  torque,  and  the 
.load  reaction  torque,  (^. 

A  PQ  -  Qa  “  Ql  ft  (16) 

.'ubstituie  equation  (5)  into  equation  (1)  and  thereby  eliminate  from 
equation  (l)  and  equate  this  result  to  to  obtain  the  following  equation. 

^  (hOfh),  Np)  -  Q]^)  rpn/sec  (17) 

Jp 
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Equatiojis  (M)  and  (17)  are  a  syatem  of  two  coupled,  nonlinear,  first- 
order,  ordinary  differential  equations,  and  hence  require  two  initial 
conditions:  Np  (0)  and  (0)  to  uniquely  detemine  the  solution. 
Consequently,  the  tine-dependent  behavior  of  a  gas  turbine  engine  is 
deducible  from  the  solution  of  the  following  initial -value-problom  (Ivp), 

Np  ”  (f  (h(Nh),  Np)  -  Q]_)  rpn/sec 
Jp 

^  JpNpNp  +  (JhNh  Jig'  (Nh))  Nh  * 

19.7  y  (Wf(t)  -  C(Qi,  Np))  lb  ft  (18) 

Np  (0)  »  Npo  rpm 

Nh  (0)  ■  Nho 

The  derivation  cf  the  initial  value  problem  (equation  (18))  is  based 
upon  certain  sliqplifying  assuoiptions  which  offer  a  ooiqiroDlse  between 
reality  and  mathematical  tractability.  It  has  been  asstused  that:  (1)  there 
Is  a  one-to-one  correspondence  between  fuel  rate  and  gas  generator  speed 
in  the  steady-state;  (2)  the  definition  of  the  specific  fuel  consumption 
can  be  extended  to  seive  as  a  meaningful  statement  of  orerell  system 
efficiency  for  both  static  and  (^rneutic  machine  operation;  (3)  the  steady- 
state  operating  characteristics  of  the  power  turbine  ean  be  used  as  a 
solution  to  the  problem  of  determining  the  aeoelerating  torque  on  the 
turbine  by  formally  replacing  with  h  (equation  (5));  (4)  the  two 
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cOBforeaaor-turblDft  pairs  chang«  speeds  in  pbase  during  a  transient;  and 
(5)  a  quasl-statio  fuel  rate  can  be  defined  to  serve  as  a  measure  of  the 
•xMrgy  xkeeded  to  sustain  the  engine  at  each  of  the  unsettled  states  that 
it  passes  throoigh  during  a  transient. 
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THE  n‘ET,  CONTROLLER 

The  fuel  controller  on  a  gaa  turbine  engine  is  a  nultl-ir.out,  single- 
outpjt  device.  Controller  response  Is  characterietioally  rapid  with  respect 
to  the  response  time  of  the  engine.  In  fnct,  the  dynamic  behavior  of  the 
fuel  oortroller  can  be  euccessfully  described  by  assiudng  that  the  equations 
governing  its  operation  do  not  expllcitely  depend  upon  time.  Consequently, 
the  tine -dependent  nature  of  the  fuel  controller  must  be  a  manifestation 
of  the  transient  behavior  of  one  or  more  of  its  several  inputs. 

No  atten5)t  will  be  made  to  analyze  the  configuration  of  indlTidual 
constituent  elements  in  a  fuel  control  system.  However,  it  is  possible  to 
construct  a  mathematical  representation  of  such  devices  from  an  analysis 
of  transient  data  acquired  by  performing  a  series  of  tests  designed  specifically 
for  this  purpose.  These  tests  are  based  on  the  assuiqption  that  the  fuel 
controller  can  be  regarded  as  a  static  nonlinearity  of  the  gas  turbine  system. 

Let  Wfn  and  be  the  maximum  and  idle  (minimum)  fuel  rates,  respectively. 
Construct  the  function,  defined  by  the  following  equation: 


W  -  W  >  *fm  -  *fl 
Vk  "fk-1  n 


Ib/hr 


h"2,,.,.,n 


(19) 


irtiere  ■  0  and  At  each  of  the  n  different  steady-state  fuel 

rates  defined  by  equations  (19),  conduct  the  following  tests  on  the  fuel 
control Isr. 
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At  t  ■  07  Wf  *  and  Njj  ■  h”^  wher®  is  the  Inverse  of  the 

function  b  defined  Implicitely  by  equation  (5).  Let  L  be  throttle  position 
and  thus  and  Ljj  correspond,  respectively  ,  to  Wfjj,  Wfi  and 

Introduce  a  step-cliange  in  L  from  to  Ljj  +  ,1  -Jsj  (m  ig  an  integer). 

ID 

Record  the  final  gas  generator  speed,  N^,  and  the  full  transient  response 
of  the  fuel  rate.  Reset  the  throttle  to  and  perform  a  test  siinilar  to 
the  previous  one  except  that  in  this  instance,  the  throttle  is  to  be 
advanced  to  the  position  Ljj  +  ^  (L^  “  Exhaust  the  m  possible  tests 

in  this  aeqxience  (Ljt  to  ^  (I®  -  L^);  J  ■  1,  .  .  .  ,  m)  and  reset  the 
throttle  to  Introduce  a  throttle  step-change  from  the  setting  Lj^  to 

Z.  hki)  .  Record  the  final,  gas  generator  speed  and  the  full 

XQ 

transient  response  of  the  fuel  rate.  Again,  exhaust  the  m  possible  tests 
in  tills  sequence  (Li£  to  L)(  -  ^  (Ljj  -  Li);  J  ■  1,  .  .  .  ,  m) .  Set  the 
throttle  to  achieve  a  fuel  rate,  ^  and  denote  this  throttle  setting 

by  I  ic  +  1  and  repeat  the  2in  tests  conducted  previously.  Continue  in  this 
maonsr  until  all  n  possible  fuel  rates  k  ■  1,  ,  .  ,  ,  n)  have  been 

explored. 

The  integer,  n,  should  be  greater  than  or  equal  to  three  (3).  The 
integer,  ra,  however,  is  evidently  dependent  on  the  difference  Ljj  - 
(Lie  •  I-i)  positive  (negative)  throttle  step-changes.  Near  (W^i) 
the  difference  Lq  -  L]c  (Li(  -  Li^  is  small  and  hence  the  number  of  tests, 
n,  required  in  this  interval  are  few  (perhaps  one).  Conversely,  large 
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diffCTti.uus  between  (Lj^)  and  L){  (Li)  for  positive  (negative)  throttle 
etep^changee  may  require  m  to  be  three  or  irore  to  enaure  that  the  transient 
records  evince  a  smooth  pattern. 

The  fuel  controller  test  records  for  the  positive  and  negative  throttle 
step  disturbances  are  annlyeed  over  the  time  interval  (0,  c)  irtiere  <  is 
some  positive  nunber  and  zero  is  the  tine  at  which  the  transient  wes 
initiated.  Let  be  the  majclnium  time  constant  of  the  fuel  controller. 

Then  ^  is  bovmded  below  by  and  above  by  the  requiranent  that  at  t  ■  <r 
gas  generator  speed,  Njj,  must  not  differ  significantly  from  its  Initial 
value.  Determine  the  maxifflum  fuel  rate  in  the  interval  (0,  c)  for  each 
transient  record.  Construct  a  graph  of  maximum  fuel-rate  vs  the  difference 
speed,  A  Nj-,,  at  each  of  the  n  steady-state  fuel  rates,  (k  ■  1,  ,  .  ,  ,  n) , 
The  quantity,  A  is  determined  for  each  record  by  the  following  prescription; 

6  Nh  “  Nh  (co  )  -  Nh  (0)  rpm  (20) 

The  time -dependent  fuel  rate  can  then  be  expressed  mathenaticelly  by  the 
following  relationship. 

Wf  -  B  {%,  A  Nh)  3b/hr  (21) 

Interpret  A  Nh  as  the  error  between  gaa  generator  epeed,  and  its 
speed  set  point,  ^^8*  Then  let  E  be  the  function  derived  from  B  which 
takes  Mbs  **h  *f* 


»f  ■  '  <"h.  "h.) 


Ib/hr 


(22) 


Deztot*  try  D  the  fvmotlon  which  relates  L  and  and  replace  in  equation 
(22) .hy  D  (L)  and  substitute  this  result  into  equation  (18)  to  obtain  the 
following  system  of  equations. 

jip  •  2*5^  (f  (h(Nh),  Np)  -  Qj^)  rp«/aeo 

i  JpNpl^  (JhHh  *  Jl  i‘  (Nh))  ^  - 

l9.7y(E  (D(L),  Hh)  -  C  (Qi,  Np))  iVft  (23) 

Hp(0)  -  Npo  rpm 

Mh(0)  ■  %o  rp® 
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CONCLIJSIOHS 

A  marlae  propulslor.  gas  turbine  engine  Is  dynamloslly  ohareoterluable 
by  a  syatein  of  two  nonlinear  differential  situations .  These  equations  are 
a  succinct  statement  of  the  phj'sical  principles  underlying  the  transient 
behavior  of  gas  turbines  in  general.  Representations  of  specific  systems 
can  be  aystematleally  deduced  from  equation  (2J)  by  ia^wslng  the  appro¬ 
priate  oathemat leal  conetrainta  on  these  equations.  These  constraints, 

In  turn,  are  developed  from  data  derived  from;  (1)  the  steady-state 
characteristics  of  the  engine,  and  (2)  the  dynamic  responses  of  the  fuel 
control  aystera  to  select  input  disturbances. 
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